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Resumen 


Álvarez-Olguín, G., 6 Escalante-Sandoval, C. A. (enero- 
febrero, 2016). Análisis de frecuencias no estacionario de 
series de lluvia anual. Tecnología y Ciencias del Agua, 7(1), 71- 
88. 


En las zonas del noroeste de México, donde existe una 
alta presión sobre el agua, las estrategias para afrontar 
periodos de escasez de este recurso podrían fallar si la 
disponibilidad se evalúa a través de valores medios de 
escurrimientos o precipitación, sin considerar el cambio 
en el comportamiento de tales variables. El objetivo de 
este trabajo fue estimar eventos de lluvia acumulada 
anual del noroeste de la república mexicana, asociados 
con diferentes periodos de retorno y escenarios futuros. 
Para siete funciones de probabilidad de valores mínimos 
se propusieron tres modelos no estacionarios válidos para 
series cuyas características estadísticas se han modificado 
a través del tiempo. Se analizaron 221 series de lluvia 
acumulada anual, de las cuales se determinó que 36 (16%) 
no son estacionarias; para 30 de éstas, los modelos no 
estacionarios representan mejor la variabilidad de los datos 
que los modelos convencionales. Se predice para el año 2044 
una disminución menor que 20% de la lluvia media anual, 
en estaciones ubicadas en la península de California, Sonora 
y Sinaloa. Los resultados sugieren que estas zonas serán más 
susceptibles a sequías en el futuro. 

Palabras clave: análisis de frecuencias, series no 
estacionarias, escasez de agua, sequías, disponibilidad de 
agua, cambio climático, noroeste de México. 


Introducción 


La predicción de eventos de lluvia mediante 
análisis de frecuencias es de suma importancia en 
diversos estudios de ingeniería, como el diseño 
de obras hidráulicas para el control de escurri- 
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mientos, la planeación adecuada de uso de suelo 
y la evaluación de la disponibilidad de agua. No 
obstante, debido a las implicaciones económicas y 
sociales, el análisis de frecuencias debe realizarse 
con alta precisión, pues la subestimación de la 
disponibilidad de agua podría implicar fallas 
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en las estrategias de asignación de agua a los 
usuarios, en particular durante periodos pro- 
longados de sequías. 

La metodología clásica de análisis y predic- 
ción de eventos de variables hidrológicas se 
basa en la teoría de valores extremos de series 
estacionarias. Bajo condiciones de estacionarei- 
dad, la distribución de la variable de interés es 
invariante en el tiempo, sin tendencias, cambios 
y periodicidades (Villarini £: Smith, 2010). No 
obstante, la validez de la estacionareidad de las 
series de tiempo debido al cambio del clima está 
siendo cuestionada por la comunidad científica 
(Milly et al., 2008; Gleick, 1989; Voss, May € 
Roeckner, 2002; Held éz Soden, 2006; Webster, 
Holland, Curry, € Chang, 2005). Según el Grupo 
Intergubernamental de Expertos sobre Cambio 
Climático (IPCC, 2007), las actividades humanas 
están provocando el calentamiento del planeta, 
por lo que se han registrado numerosos cambios 
de largo plazo en el clima, que incluyen posi- 
bles aumentos en la intensidad de los ciclones 
tropicales y las ondas de calor, así como en la 
intensidad y frecuencia de eventos extremos, co- 
mo sequías y lluvias torrenciales. Se espera que 
el cambio climático intensifique el estrés actual 
sobre los recursos hídricos debido al crecimiento 
de la población, a los cambios económico y de 
uso de la tierra, y en particular a la urbanización 
(IPCC, 2007). Por lo anterior, conocer la disponi- 
bilidad del agua es un factor clave en la creación 
de estrategias para afrontar periodos de escasez, 
no obstante, el proceso de asignación de agua a 
los usuarios podría ser ineficiente si la disponi- 
bilidad se evalúa a través de valores medios de 
escurrimientos o precipitación, sin considerar el 
cambio en el comportamiento temporal de estas 
variables. Por el contrario, si la disponibilidad 
se estima a partir de eventos de lluvia asociados 
con una determinada probabilidad y escenarios 
que consideren la disminución de la lluvia, se 
podrán crear mejores estrategias para la distri- 
bución del agua. 

En México, los escenarios climáticos cons- 
truidos hasta ahora sugieren que la precipitación 
en muchas regiones del país podría disminuir; 
por lo tanto, bajo un escenario de cambio cli- 
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mático con condiciones de sequía prolongada, 
el desarrollo sostenido de la sociedad estaría en 
peligro (Magaña-Rueda éz Gay-García, 2002); 
principalmente en zonas del noroeste del país, 
donde existe una alta presión sobre los recursos 
hídricos (Conagua, 2012). Ante este panorama, 
el objetivo de este trabajo fue estimar eventos 
de lluvia acumulada anual del noroeste de la 
república mexicana, asociados con diferentes 
periodos de retorno y escenarios futuros, a 
través de modelos probabilísticos válidos para 
series de tiempo no estacionarias. Los resultados 
obtenidos podrán tomarse como base para una 
adecuada evaluación de la disponibilidad de 
agua en la región. 


Materiales y métodos 


La zona de estudio abarcó los estados de Baja 
California Norte, Baja California Sur, Sonora, 
Sinaloa, Chihuahua y Durango, ubicados en 
el noroeste de la república mexicana (figura 
1). Estas entidades fueron seleccionadas por 
tener altas presiones sobre los recursos hídricos 
(Conagua, 2012) y zonas con alta vulnerabilidad 
ante las sequías. 

Los datos pluviométricos utilizados en este 
trabajo corresponden a los registros diarios de 
las estaciones climatológicas convencionales 
ubicadas en la zona de estudio, contenidas en 
la base de datos CLICOM (Climate Computing 
Project) del Servicio Meteorológico Nacional 
de la Comisión Nacional del Agua (Congua). 
Se consideraron estaciones con al menos 90% 
de registros completos entre el periodo de 1950 
a 2013. Mediante la prueba de Grubbs y Beck 
(1972) se identificaron valores de lluvia anual 
que se desvían de modo considerable de la dis- 
tribución de los datos, y fueron comparados con 
los registros de estaciones aledañas para corro- 
borar su magnitud. Los valores que resultaron 
claramente erróneos fueron eliminados. Para 
deducir datos faltantes se aplicó el método de 
interpolación de la distancia inversa ponderada 
(Shepard, 1968), en el cual se consideraron dos 
estaciones de apoyo y un exponente de distancia 
igual a dos. Para verificar si existe mejora en 
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Figura 1. Localización del área de estudio. 


la estimación de la media y la varianza de las 
muestras extendidas, se utilizó el criterio de 
información relativa (1), definido por Escalante- 
Sandoval y Reyes-Chávez (2002) como: 


me Var(S;,) 
Var(S ) 


2 
y 


(1) 


Donde Var(S; ) es la varianza de la varianza 
de la serie original, y Var(s;,) es la varianza de 
la varianza de la serie extendida. 

Si 1 > 1, entonces la varianza de la varianza 
de la serie extendida no excede a la original, por 
lo que la extensión de los registros es adecuada. 
Con base en lo anterior, se obtuvieron datos 
confiables de lluvia diaria de 221 estaciones. 

Para poder comparar espacialmente los cam- 
bios ocurridos en las series de lluvia, se llevó a 
cabo una regionalización de la lluvia mensual 
basada en un análisis de componentes prin- 
cipales (ACP) en modo S (Richman, 1986). Se 
seleccionaron las primeras cinco componentes 
que explican 73% de la varianza de la precipita- 
ción. Las componentes retenidas se sometieron 
a una rotación oblicua directa, con el parámetro 


de oblicuidad d = 0. En la delimitación de las 
regiones homogéneas se asignó a cada estación 
la componente principal de mayor carga en 
valor absoluto, considerando sólo aquellas con 
valores mayores a 0.4. Las regiones se definieron 
de acuerdo con el área de influencia de cada 
estación obtenida con el método de polígonos 
de Thiessen. 

El siguiente paso fue verificar la presencia 
de cambios graduales en las series (tendencias) 
o puntos de cambio abruptos en la media y/o 
la varianza de la distribución de la variable de 
interés (Villarini et al., 2009). Se aplicó la prueba 
de Pettitt (Pettitt, 1979), que permite detectar 
cambios significativos en la media (mediana) 
de las series cuando el punto de cambio es 
desconocido. Las tendencias en las series se 
analizaron con el método de descomposición 
de modos empíricos (DME) desarrollado por 
Huang et al. (1998), el cual es un método intui- 
tivo, directo y adaptativo a procesos no lineales 
y no estacionarios, cuya efectividad para el 
análisis de tendencias y señales de periodo largo 
se ha demostrado en estudios recientes (Oian, 
Wu, Fu, € Zhou, 2010; Carmona, 2010; Wu, 
Huang, Long, €: Peng, 2007). Para determinar 
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la significancia de las tendencias, los residuales 
obtenidos por DME se analizaron con la prueba 
de Mann-Kendall (Kendall, 1975). 

De los fenómenos de gran escala que influ- 
yen en la variabilidad climática en México, se se- 
leccionó la oscilación decadal del Pacífico (PDO, 
por sus siglas en inglés) para ser incluida en el 
análisis de frecuencias de la lluvia, esto debido 
a la asociación existente entre episodios persis- 
tentes de valores PDO negativos y condiciones 
de sequías severas y prolongadas a través del 
occidente y centro de América del Norte (Cook, 
Woodhouse, Eakin, Meko, éz Stahle, 2004). Por 
tal razón, se analizó la correlación entre las se- 
ries anuales regionales de lluvia acumulada y 
el índice PDO, a través del coeficiente de corre- 
lación p (rho), propuesto por Spearman (1904), 
el cual es una prueba no paramétrica que mide 
la asociación entre dos variables discretas y está 
definido por la ecuación (2): 


2 
pato E 0) 
N(N? -1) 

Donde los valores de p varían en el intervalo 
[-1, 1]; d es la diferencia entre los valores corres- 
pondientes de x-y; N, el número de parejas, para 
N >20, sigue una aproximación a la distribución 
t de Student (3): 


(3) 


En el análisis de frecuencias, para conside- 
rar la falta de estacionareidad de la series se 
pueden asociar los parámetros de las distribu- 
ciones con más de una covariable. La inclusión 
de covariables ha permitido incorporar 
tendencias o índices de circulación atmosférica 
de gran escala como ENSO (El Niño-Southern 
Oscillation), PDO (Pacific Decadal Oscillation 
o NAO (North Atlantic Oscillation), en la 
modelación de datos de precipitación extrema 
(Katz, Parlange, £: Naveau, 2002; El Adlouni, 
Ouarda, Zhang, Roy, €: Bobée, 2007; Villarini, 
Smith, £ Napolitano, 2010; Vovoras €: Tsokos, 
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2009; Ouarda éz El Adlouni, 2011), temperaturas 
máximas (Villarini et al., 2010), gastos máximos 
(Villarini et al., 2009; Towler et al., 2010; Vogel, 
Yaindl, €z Walter, 2011), nivel extremo del mar 
(Coles, 2001) e intensidades de huracanes (Mes- 
tre £: Hallegatte, 2009). 

No existe una teoría general establecida 
para procesos no estacionarios, por lo que es 
usual utilizar como base modelos de valores 
extremos estándares y extenderlos (Coles, 2001). 
Por ejemplo, de la distribución log normal (LN3) 
con parámetros x, u y O (ubicación, escala y 
forma, respectivamente), se puede obtener un 
modelo adecuado para estimar una variable x, 
en función del tiempo, de la siguiente forma: 


x, - LN3(x,(4),u(),0(1)) (4) 


Donde cada parámetro x,(+), u(t), o(t) tienen 
una expresión en términos del tiempo. El pará- 
metro u(t) se puede expresar con la ecuación (5): 


u(t) = 1, +u,t (5) 


Donde t corresponde a la tasa anual de cam- 
bio en el tiempo. De esta forma, variaciones a 
través del tiempo en los procesos observados 
se modelan como una tendencia lineal en el 
parámetro de ubicación de la distribución. 

La estimación de los parámetros de un 
modelo en función del tiempo, o de otra cova- 
riable, se puede realizar mediante el método de 
máxima verosimilitud (Katz et al., 2002), consi- 
derado como el más eficiente, pues proporciona 
la menor varianza muestral de los parámetros 
estimados y, por ende, de los eventos estimados, 
en comparación con otros métodos; además, 
una ventaja es su adaptabilidad a cambios en 
la estructura del modelo. 

Para un modelo no estacionario propuesto 
de la función LN3, la función de verosimilitud 
es la siguiente: 
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Donde x,(t), u(t) y o(t) son los parámetros de 
ubicación, escala y forma, respectivamente, en 
función del tiempo o de otra covariable. 

En el análisis de lluvia acumulada anual del 
noroeste de México es de interés conocer cuáles 
serán los escenarios más adversos de escasez 
de agua en el futuro. Por lo tanto, la predicción 
de eventos de lluvia se debe realizar con fun- 
ciones de probabilidad para valores mínimos. 
Sin embargo, no existen modelos probabilísticos 
no estacionarios para funciones de distribución 
de valores mínimos. Por tal razón, se propu- 
sieron los modelos descritos en el cuadro 1. 
Los modelos M, corresponden al análisis de 
frecuencias convencional, mientras que M, M, 
y M, corresponden al caso no estacionario y 
sus parámetros se expresaron en términos del 
tiempo y del índice PDO. La serie promedio 
anual del índice PDO se obtuvo a partir de los 
valores mensuales del periodo de 1950 a 2013, 
descargados de la página del National Climatic 
Data Center de los Estados Unidos de América 
(NCDC, 2013). 

Para cada estación, los estimadores de los 
parámetros de las funciones se obtuvieron 
tras maximizar las funciones logarítmicas de 
verosimilitud respectivas. Los modelos que 
mejor describen la variabilidad de los datos se 
seleccionaron mediante el Criterio de Informa- 
ción de Akaike (1974), en el cual se utilizó la 
expresión (7). El mejor modelo es el que tuvo 
el menor valor AIC. En el caso de tener valores 
similares de AIC, se seleccionó el modelo con 
menos parámetros: 


AIC =-2(1) + 2K (7) 


Donde l es el máximo valor de la función de 
verosimilitud y K es el número de parámetros 
estimados. 

Se utilizó el método descrito por Coles (2001) 
para comparar la validez de un modelo M, 
contra otro M,, tal que M, CM, en el cual se 
utiliza la medida de discordancia definida por 
la expresión (8). 


D=241, (M,)-1,(M,)) (8) 


Donde I(M,) es el máximo valor de la función 
logarítmica de verosimilitud del modelo M.. 

El estadístico D se distribuye de acuerdo con 
la distribución chi-cuadrada (x?). El parámetro 
v es la diferencia entre el número de paráme- 
tros de los modelos M, y M,. Valores de D más 
grandes que las cantidades de la distribución x?, 
para un nivel de confianza particular se consi- 
deran significantes; entonces se rechaza M, en 
favor del modelo M.. 

Una vez seleccionado el mejor modelo fue 
necesario verificar que se ajustara bien a los 
datos observados. Debido a que en una serie 
no estacionaria los datos no se distribuyen de 
manera idéntica, éstos deben ser transformados 
(Coles, 2001; Katz et al., 2002). La transformación 
de los datos consistió en la normalización de los 
mismos, proceso al que Krzysztofowicz (1997) 
denominó como transformación normal de 
cuantiles. La normalización se realizó mediante 
la expresión (9), en la que una variable aleatoria 
X, con una determinada distribución de proba- 
bilidad F, se transforma en una variable Y, con 
una distribución normal estándar G: 


Y =G” (F(Xx)) (9) 


Donde G” es la función inversa de la distri- 
bución normal estándar definida por la ecuación 
(10): 


G(2)===S+? dz (10) 


De esta forma, tras ordenar los valores estan- 
darizados y asociarlos con los correspondientes 
de la distribución empírica de Weibull, se obtu- 
vieron los siguientes pares de datos: 


da mai 
(%,/0 (1) jo... (11) 


Con los pares de datos anteriores se elabo- 
raron las gráficas de cuantil-cuantil (Q-Q) y 
“worm plot” (Buuren éz Fredriks, 2001). En una 
gráfica “worm plot”, el eje vertical es la dife- 
rencia entre los valores empíricos y teóricos, y 
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contiene el intervalo de confianza a 95% por la 
ecuación (12): 


51.96 (2) ' (p(1-p)/m) (12) 


Donde Á(z) es la función de densidad normal; 
z es un evento asociado con una probabilidad, y 
n es el tamaño de la muestra. 

Finalmente, se estimaron eventos de diseño 
X asociados con diferentes periodos de retor- 
no. Para la función LN3, con el modelo M,, los 
eventos de diseño se estimaron al resolver la 
ecuación (13): 


X=3,+exp(F" (p/6(4),5)) (13) 
ñ(+) =6,+ñ,t,+6,PDO (14) 


2-ñ(1) 


perio a 


0y2x1 


Donde x, es el estimador del parámetro de 
ubicación; 6, el estimador del parámetro de 
forma; ú(t), el parámetro de escala en función 
de los estimadores '',, ", y ",; £, el valor del ín- 
dice del tiempo para un escenario determinado 
(t,= 65 para el escenario del año 2014 y t, = 95 
para el escenario del año 2044); PDO, un valor 
promedio del índice PDO, y X es el evento de 
diseño asociado con un periodo de retorno T y 
un escenario futuro t,. 

El análisis de datos abarcó el periodo com- 
prendido entre los años 1950 y 2013, por lo cual, 
al primer año (1950) se le asignó el valor t,= 1, 
mientras que al último (2013), le correspon- 
dió un valor t,= 64. De esta forma, para años 
posteriores, los valores del índice de tiempo 
aumentaron una unidad de forma consecutiva. 

Debido a la asociación de valores PDO ne- 
gativos y periodos con déficit de lluvia, el valor 
PDO se estimó como la media aritmética de los 
valores PDO correspondientes al lapso de ano- 
malías con mayor déficit de lluvia acumulada 
anual analizada. 


Resultados y discusión 


La distribución espacial de la lluvia acumulada 
anual en la zona de estudio se encuentra en la fi- 
gura 2. Los valores más bajos se presentan en la 
península de California, con menos de 200 mm 
al año, siendo la zona con menor precipitación 
de todo el país; en contraste, entre los límites de 
Sonora, Sinaloa, Chihuahua y Durango se tienen 
valores de más 1 000 mm. 

El análisis de componentes principales 
arrojó cinco componentes que explican 73% 
de la varianza de la precipitación. Tras asignar 
la componente de mayor carga al polígono 
de Thiessen de cada estación, se pudieron 
identificar espacialmente grupos de polígonos 
asociados con la misma componente; cada gru- 
po se consideró como una región homogénea; 
de esta forma, se generaron las cinco regiones 
mostradas en la figura 2. 

De acuerdo con el cuadro 2, 19% de las series 
analizadas es no estacionaria, de las cuales, 47% 
se ubica en la región 1, que abarca principal- 
mente los estados de Sonora y Sinaloa. Del total 
de series que presentan puntos de cambio des- 
cendentes significativos, 73% corresponde a la 
región 1; en esta región, además, se localizan las 
tres series que presentan puntos de cambio des- 
cendentes y tendencias decrecientes, así como 
dos de las tres series con tendencia decreciente. 

Según el análisis de correlación de Spear- 
man, existe una correlación de 0.43, estadísti- 
camente significativa a cualquier nivel, entre 
el índice PDO y la serie de lluvia acumulada 
de la región 3; mientras que con las regiones 4 
y 5, la correlación es de 0.2, significativa a un 
nivel de 0.1. Para las tres regiones anteriores, la 
correlación es positiva; por lo tanto, al estar el 
PDO en fase negativa de manera persistente, se 
tendrán periodos con déficit de lluvia, tal como 
ocurrió en la región 3 entre 1950 y 1975 (figura 
3a). La asociación entre episodios prolongados 
de valores PDO negativos y condiciones de 
sequías severas y prolongadas a través del 
occidente y centro de Norteamérica fue men- 
cionada por Cook et al. (2004). Entre 1977 y 1997 
se presentó un régimen PDO positivo y apa- 
rentemente a partir de 1998 se tiende hacia una 
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Figura 2. Regiones homogéneas y distribución de la lluvia acumulada anual. 


Cuadro 2. Número de series de lluvia acumulada anual, según su condición determinada por las pruebas de Pettitt 


y Mann-Kendall, a un nivel de significancia de 0.05. 


Región Total E NE A D TC TD ATC De€TD 
1 72 55 17 2 8 0 2 2 3 
2 33 27 6 1 0 5 1 1 0 
3 43 38 5 Z 3 0 0 0 0 
4 34 28 6 0 0 2 0 4 0 
5 39 37 2 1 0 0 0 1 0 
Total 221 185 36 6 11 5 3 8 3 


E: estacionaria; NE: no estacionaria; A: cambio ascendente; D: cambio descendente; TC: tendencia creciente; TD: tendencia decreciente; ATC: 
cambio ascendente y tendencia creciente; Dé*TD: cambio descendente y tendencia decreciente. 


fase negativa, por lo que es muy probable que 
en las próximas décadas cambie de polaridad 
(Mantua, Hare, Zhang, Wallace $z Francis, 1997; 
MacDonald € Case, 2005). Debido a lo anterior, 
es posible que la pérdida de la estacionareidad 
de las series por el incremento de la lluvia se 
deba en gran medida a la influencia de la fase 
positiva del PDO, por lo tanto, la persistencia 
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de valores PDO negativos que pudieran causar 
sequías prolongadas se consideró en el análisis 
de predicción de la lluvia. 

Del análisis de frecuencias realizado a las 221 
series de lluvia acumulada anual, las pruebas de 
bondad de ajuste indicaron que para 84 de éstas, 
los modelos no estacionarios estadísticamente 
representan mejor la variabilidad de los datos 
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(cuadro 3). En el caso de las series no estaciona- 
rias, en 83% de los casos los modelos no estacio- 
narios se ajustó mejor a los datos; para las series 
no estacionarias restantes, aunque el estadístico 
AIC indicó que los modelos no estacionarios fue- 
ron mejores, la medida de discordancia (D) fue 
inferior a 3.84 (valor de yx? a un nivel a. = 0.05), 
por lo que no se justificó el uso de un modelo M,, 
M, o M, en lugar de un M,. 

La mayoría (60%) de los modelos no estacio- 
narios seleccionados es M, o M, de la función 
LN3. Los resultados obtenidos son consistentes 
con las pruebas de estacionareidad y de corre- 
lación realizadas; por ejemplo, para la región 3, 
que tiene la correlación más alta con PDO, se 
seleccionaron modelos M, en 20 casos. Además, 
para la región 1, donde se localiza 47% de las 
estaciones con series no estacionarias, se selec- 
cionaron 22 modelos M,. Es importante señalar 
que el modelo M, de la distribución GW se selec- 
cionó para 61 de las series estacionarias; por lo 
tanto, las distribuciones mezcladas representan 
mejor la variabilidad de los datos cuando éstos 
son estacionarios. 

Como ejemplo del procedimiento para la 
selección del mejor modelo, en los cuadros 4 
y 5 se muestran los resultados de las pruebas 
de bondad de ajuste y los estimadores de los 
parámetros de los modelos ajustados a las series 
de las estaciones 2060 (“Santa Cruz”, Ensenada, 
Baja California) y 25082 (“San Blas”, El Fuerte, 
Sinaloa), respectivamente. Cabe mencionar que 
la serie de la estación 2060, según la prueba de 
Pettitt, presentó en el año 2002 un punto de 


quiebre significativo a un nivel de 0.05, con 
un cambio descendente en la media de la serie 
(figura 3b). Por otro lado, para la estación 25082, 
el residual obtenido por el método DME (figura 
3c) indicó que existe una tendencia monótona 
decreciente, estadísticamente significativa la 
un nivel de 0.05), según la prueba de Mann- 
Kendall. 

Para la estación 2060, el modelo con el menor 
valor AIC fue el M, de la función LN3, en el cual 
el parámetro de escala está en función del índice 
PDO; además, el valor del estadístico D es supe- 
rior a 3.84, lo cual indica que es más adecuado 
y explica mejor la variabilidad de los datos 
que el modelo M,,. Para la serie de lluvia de la 
estación 25082, el modelo M, de la función GA3, 
que incluyó una tendencia en el parámetro de 
ubicación, representa mejor la variabilidad de 
los datos que los modelos restantes; el valor de 
AIC del modelo M, fue el menor, con un valor D 
igual a 5.0, estadísticamente significativo. 

La figura 4 corresponde a las gráficas “worm 
plot” y de cuantil-cuantil (Q-Q) de las estaciones 
2060 y 25082, para los modelos mencionados. Se 
observa para ambos casos que los datos de las 
gráficas caen dentro de los límites de confianza y 
están cerca de la diagonal unitaria; por tal razón, 
se determina que los modelos seleccionados se 
ajustan de forma adecuada a los datos. 

Para la estación 2060, el valor del índice PDO 
utilizado para estimar los eventos de diseño fue 
de -0.86, promedio del periodo comprendido en- 
tre 2007 y 2013, caracterizado por la persistencia 
de valores negativos y que abarcó el periodo 


Cuadro 3. Número de modelos seleccionados para la predicción de eventos de lluvia acumulada anual por región. 


Log normal 3p Gamma 3p Weibull ne 
Weibull 
Región | Total M, M, M, M, a M, M, M, M, M, 
1 12 12 7 27 
2 33 2 6 6 
3 43 0 2 5 
4 34 a 2 9 
5 39 1 6 14 
Total 221 18 23 61 
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Figura 3. Variación temporal de series anuales del noroeste de México: a) lluvia media de la región 3 e índice PDO); b) estación 
2060, “Santa Cruz”, Ensenada, Baja California; c) estación 25082, “San Blas”, El Fuerte, Sinaloa. 


con mayor déficit de lluvia (figura 3b). Por lo Los valores de los eventos de diseño asocia- 
anterior, para el modelo M, de la función LN3 dos con diferentes probabilidades de no exce- 
(cuadro 4), los estimadores de los parámetros dencia, p, se estimaron al resolver la siguiente 
fueron los siguientes: expresión: 


X 0.05 +expfLN3” (p[5.40,0.43)) 
ñi(+) = 4, +1,PDO =5.65+(0.29)(-0.86) 
Donde LN3” es la función de distribución de 


a o probabilidad log normal inversa. 
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Figura 4. Gráficas “worm plot” (izquierda) y Q-Q (derecha) para el análisis visual del ajuste de los modelos seleccionados: a) 
estación 2060, “Santa Cruz”, Ensenada, Baja California; b) estación 25082, “San Blas”, El Fuerte, Sinaloa. 


Por otro lado, para la estación 25082 (cuadro Las expresiones utilizadas para estimar los 
5), los parámetros estimados para el modelo M, eventos de diseño fueron las siguientes: 
de la función GA3 fueron los siguientes: 


Para el escenario del año 2014: 
Para el escenario del año 2014: 


a X =-408.82 + GA3” (p]26.43, 33.17) 
2,(0 =2, +4, f,=-240.47 + (2.59)(65) 
Si Para el escenario del año 2044: 


Para el escenario del año 2044: A 
X = -486.52 + GA37 (p| 26.43, 33.17) 


2, (1) =2, +2, t, =-240.47 + (2.59)(95) A 
— 486.52, 38.17, $ = 26.43. onde GA3” es la función de distribución de 


probabilidad gamma inversa. 
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Cuadro 4. Resultados de las pruebas de bondad de ajuste para la estación 2060 (“Santa Cruz”, Ensenada, Baja California). 


Estimadores de los parámetros 


Estación | Función | AIC D 


iy Ñ, la %a ó 
2060 | LN3(M,) | 8153 | - 5.61 - - 0.00 [050 | - - - - - | - 
2060 | LN3(M) | 817.3 | 0.0 | 5.63 0.00 - 0.00 [0.50 | - - - - - | - 
2060 | LN3(M,) | 799.8 | 17.5 | 5.65 0.29 - 0.05 [043 | - - - - - | - 
2060 | LN3(M,) | 7999 | 19 | 5.79 0.00 | 0.31 | 0.02 [0.43 | - - - - 0 E 
> di e a p 
2060 | GA3(M,) | 815.6 | - 64.50 - -  |10256|2.37 | - - - - - | - 
2060 | GA3(M) | 816.0 | 16 | 84.02 | -0.42 - | 103.77 (229 | - . , a 2 |. 
2060 | GA3(M2) | 810.3 | 7.3 | -1971 | 56.74 | - | 5417 |621 | - - - - | - 
2060 | GA3(M3) | 811.0 | 13 | 5.87 -0.81 | 59.58 | 52.71 [640 | - - E a Sl 
Y % A p a 
2060 | W3(M) | 816.5 | - 77.14 - - | 33228 |148| - - - - 0 E 
2060 | W3(M) | 8120 | 65 | 133.29 | -0.91 -  |318.37|122| - - ó a 0] 
2060 | W3(M,) | 8183 | 02 | 7200 | -102 - | 32685 |149| - - - - | - 


2060 W3 (M,) | 808.1 | 12.2 168.09 -1.40 16.57 | 323.54 | 0.98 - - - - - - 


0 l 2 
2060 G(M,) 870.0 = 397.17 = - 207.06 - = = - - = - 
2060 G(M,) 870.9 | 1.1 334.58 1.93 - 203.70 - - - - - - - 
2060 G (M,) 846.2 | 25.8 | 407.61 146.91 - 162.28 - - - - - - - 


2060 G(M,) 847.5 | 0.7 445.44 -1,.14 | 154.50 | 161.04 5 = > = a Z 5 


d Ú, Ú, Br á, d, D, d, P, a |? 
2060 | ww(M,) | 822.2 772 : O O O O o O 
“9 2060 | ww(m) | 8071 191| 1332 | -09 | - | 3184 | 16 [1345 |-09| - |3183|07 | 07 
E 2060 | ww(m) si66| 96 | 725 | 105 | - | 3271122 | 83.4 |-03 | - |352|05|08 
5 2060 | WW(M) 8297 |-91 | 1675 | -24 | 173 | 3234 | 10 |167.5 | -14 | 163 [3234 | 10 | 0.6 
E R, £, 3 á, ña a. a/a A 
5 2060 | GG(M) | 8394 | - | 528 k O a A esa E a 
- 2060 | cc(m) ss97|37 | 2040 | mi | - | 1523 | - |oss3|02 | - |895|- |03 
; 2060 | cc) |si6z 272 | 5780 | 2077 | - | 1065 | - |30m6|s1] - | 776 | - [oz 
E 2060 | cc(m) sizo|23 | 62.1 | -26 | 2800 | 1230 | - |3185|-07 | 673 | 660 | - |03 
a R IÓ a, a lolo l|salal> 
E 2060 | cw(m) |825| - | 3972 E 0 O 0 A O e 25 (154 Fo 
a 2060 | cw(m) sie7|7s | 2843 | 149 | - | 1201 | - [1307 -09| - [|3042|13 01 
E 2060 | cw(m) |s201| 64 | 5743 | 2054 | - | 1mss | - |659|-856| - |2015| 18 | 02 
E 2060 | cw(m) |sis3 58 | 454 | -11 | 1545 1545 | - [168.1 |-14 166 | 3235 | 10 | 01 
S AIC: criterio de información de Akaike (1974) (adimensional); D: medida de discordancia (adimensional). 
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Cuadro 5. Resultados de las pruebas de bondad de ajuste para la estación 25082 (“San Blas”, El Fuerte, Sinaloa). 


Estimadores de los parámetros 


Estación | Función AIC D 


hy 6, | 6, % 6 
25082 | LN3(M,) | 852.1 6.26 - - 0.00 | 034 | - - - - - | - 
25082 | LN3(M,) | 848.9 | 51 | 643 |-0.01| - | -0.01 | 033 |  - - - - - | - 
25082 | LN3(M,) | 854.0 | 0.0 | 626 | 0.00 | - 0.00 | 034 | - - - - - | - 


25082 LN3(M,) | 850.5 | 0.4 6.44 | -0.01 | 0.03 0.02 0.33 = = = > = = 


De E ó B 
25082 | GA3(M,) | 848.7 -241.20 | - - | 39.94 | 1986 |  - - - - - | - 
25082 | GA3(M) | 8458 | 5.0 |-240.47 | -259| - | 33.17 | 26.43 |  - z E : eN ME 
25082 | GA3(M,) | 850.7 | 0.0 |-241.20 | -0.88 | - | 39.95 | 19.85 |  - - - - - | - 
25082 | GA3(M,) | 847.2 | 0.6 | -990.42 | -2.98 | 21.56 | 17.22 | 95.36 |  - ] - - EL 
Yo A | p a 
25082 | W3(M,) | 848.4 125.68 | - - | 605.57 | 257 | - - - - - | - 
25082 | W3(M) | 8502 | 02 | 19718 |-110| - [607.81 | 233 | - : E , 0 E 
25082 | W3(M,) | 847.4 | 3.0 | 89.13 |-110.0| - |60324 | 272 |  - - - - - | - 


25082 W3(M,) | 849.3 | 3.0 9.01 1.16 | -133.4 | 600.10 | 2.96 5 a ss = a = 


D O, O, Gá 
25082 G(M,) 866.0 643.63 - - 191.84 - - - - - - - 
25082 G(M,) 857.2 | 10.8 | 768.42 | -4.01 - 172.22 - - - - - - - 
25082 G(M,) 865.5 | 2.5 | 650.44 | 44.64 - 185.02 - - - - - - - 


25082 G(M,) 854.2 | 4.9 | 776.41 | -4.03 | 56.49 | 162.50 - - - - - - - 


% d, d, pr Ó, 0, b, d, A 

25082 WW(M,) | 852.9 124.9 - - 622.4 8.5 139.5 - - as | 22 | 02 
25082 WW(M,) | 856.3 | 0.6 194.3 -0.3 - 621.7 73 | 207.2 | -1.1 - 593.2 | 2.0 | 0.2 ES 
25082 WW(M,) | 855.4 | 15 89.1 - -109.8 | 603.2 2 89.1 - -110.0 | 603.2 | 2.7 | 0.2 E 
25082 WW(M,) | 858.9 | 1.3 6.5 2.5 |-1278| 593.9 3.0 11.5 0.4 |-1348| 5972 | 2.8 | 0.6 S 
E 7 ña lala ja p 2 
25082 GG(M,) | 854.7 810.4 - - 171.4 - 561.7 - - 114.9 - 0.2 E 
25082 GG(M,) | 863.2 | -4.5 7684 | -3.9 - 172.2 - 768.4 | -4.0 - 1722 | - 0.1 E 
25082 GG(M,) 853.3 | 5.4 886.1 = 187.2 | 65.8 - 578.2 - -4.4 | 126.7 | - 0.1 E 
25082 GG(M,) 847.6 | 195 | 979.8 -4.1 | 1163 | 63.0 - 675.8 | -3.2 16.8 | 111.1 - 0.2 á 
R, RR 0, 0, o | 0, E = 
25082 GW(M,) | 850.9 625.2 - - OJD - 139.4 - - NO | 22 | 02 E 
25082 GW(M) | 852.2 | 2.7 632.2 75 - 57.2 - 106.5 | -109.0 - 594.5 | 2.4 | 0.2 » 
25082 GW(M,) | 8522 | 2.7 632.2 - YES) 02 - 106.5 - -109.0 | 594.5 | 2.4 | 0.2 > 
25082 GW(M,) | 855.8 | 0.4 630.7 | -8.0 | -15 49.2 - 109.9 0.8 -95.6 | 604.8 | 2.6 | 0.0 E 
AIC: criterio de Información de Akaike (1974) (adimensional); D: medida de discordancia (adimensional). S 
3 
a 
E 
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Figura 5. Eventos de diseño de la lluvia acumulada anual: a) estación 2060, “Santa Cruz”, Ensenada, Baja California; b) estación 
25082, “San Blas”, El Fuerte, Sinaloa. 


La representación gráfica de los eventos de 
diseño estimados para las estaciones 2060 y 
25082, a partir de las expresiones anteriores, se 
muestran en la figura 5. 

Para la estación 2060, los eventos estimados 
con una probabilidad de ocurrencia de 50% (pe- 
riodo de retorno de dos años) fueron de 271.8 
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mm, con el modelo estacionario, y 266.8 mm 
con el modelo no estacionario. Con el enfoque 
tradicional, los eventos están sobreestimados, 
mientras que con el análisis no estacionario se 
predijeron valores de lluvia para las condiciones 
históricas más adversas de las sequías originadas 
por la persistencia de valores PDO negativos. 
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Al considerar la presencia de una tendencia 
decreciente de la lluvia en la estación 25082 fue 
posible simular el decremento en la magnitud 
de los eventos de diseño para escenarios futu- 
ros. Por ejemplo, para un periodo de retorno 
de dos años se predice una disminución de la 
lluvia de 77 mm en el año 2044 con respecto a 
la lluvia del escenario de 2014. Debido a que las 
pruebas estadísticas aplicadas a la serie de esta 
estación mostraron una tendencia decreciente 
significativa, la predicción del decremento en 
la lluvia para el escenario 2044 es consistente 
con el comportamiento de la serie. 

De las demás estaciones cuyos modelos 
seleccionados (M, o M,) están en función del 
tiempo, se estimó el incremento (decremento) 
de la lluvia que se tendrá para el año 2044 
con respecto al año 2014, para un periodo de 
retorno de dos años. Este periodo de retorno 


corresponde al de la lluvia mediana (lluvia aso- 
ciada con 50% de probabilidad de ocurrencia). 

En la figura 6 se muestra la variación espa- 
cial del incremento (decremento) de la lluvia, 
como porcentaje de la lluvia media anual; para 
las estaciones que no presentan cambios, los 
eventos de diseño se obtuvieron con mode- 
los M, o modelos M,, que no dependen del 
tiempo. En 29 estaciones de las regiones 1, 2, 
3 y 5, ubicadas en la península de California, 
Sonora y Sinaloa, se predice una disminución 
menor a 20% de lluvia media anual (cuadro 6). 
Cabe mencionar que en esta zona se tienen los 
valores de lluvia media anual más bajos del 
país, por lo cual estos resultados sugieren que 
será más susceptible a sequías en el futuro. En 
contraste, en ocho estaciones de la región 4 se 
esperan aumentos en la lluvia media de hasta 
30%. 


115 W 110 W 


30 N 


Límite estatal 


[7] Regiones homogéneas 
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Figura 6. Incremento (decremento) de la lluvia entre el escenario de 2014 y 2044, con respecto a la lluvia media anual, para un 


periodo de retorno de dos años. 
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Cuadro 6. Número de estaciones según los rangos del porcentaje de incremento (decremento) de la lluvia entre el escenario 


de 2014 y 2044, con respecto a la lluvia media anual, para un periodo de retorno de dos años. 


Incremento (decremento) Región 

o 1 2 3 4 5 Total 

20 a-10 5 1 2 0 0 8 

-10a-0 14 0 5 0 2 21 

0 48 27 33 26 35 169 

0-10 2 1 3 2 1 9 

10-20 1 1 0 3 1 6 

20-30 2 2 0 1 0 5 

30-40 0 1 0 pd 0 3 

Total 72 33 43 34 39 221 

Conclusiones Con la inclusión de una tendencia en el 


En 83% de las series no estacionarias analizadas 
los modelos seleccionados para la estimación de 
eventos de diseño de lluvia acumulada en el no- 
roeste de México incluyó una tendencia lineal y 
al índice PDO como covariable; por lo tanto, se 
determina que ante la falta de estacionareidad 
en las series de lluvia, los modelos no estaciona- 
rios propuestos explican mejor la variabilidad 
de los datos que los modelos convencionales. 
Para las series no estacionarias restantes, los 
modelos no estacionarios fueron mejores; sin 
embargo, la medida de discordancia indicó 
que para un nivel de significancia de 0.05, no 
se justifica el uso de un modelo no estacionario 
en lugar de un convencional. Por tal razón, 
el cambio en el comportamiento temporal de 
la lluvia no es una condición para todas las 
estaciones. 

Los modelos no estacionarios permitieron la 
estimación de eventos de diseño asociados con 
las condiciones más adversas del índice PDO 
que han influido en los periodos históricos con 
déficit de lluvia más severos. Es posible que en 
las próximas décadas se presente de nuevo un 
periodo persistente de sequías, pues aparente- 
mente a partir de 1998 se tiende hacia una fase 
negativa del PDO) por lo tanto, la persistencia 
de valores PDO negativos deben considerarse 
en los análisis de predicción de la lluvia. 
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análisis de frecuencias se obtuvieron eventos 
de diseño asociados con diferentes periodos 
de retorno y escenarios futuros, con los cuales, 
en 29 estaciones ubicadas en la península de 
California, Sonora y Sinaloa, se predice una 
disminución para al año 2044 de menos de 20% 
de la lluvia media anual. En esta zona se tienen 
los valores de lluvia media anual más bajos del 
país, por lo cual tales resultados sugieren que 
será más susceptible a sequías en el futuro. 
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